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ABSTRACT 

Active galactic nuclei (AGN) have been observed to vary stochastically with 10 — 20% rms ampli- 
tudes over a range of optical wavelengths where the emission arises in an accretion disk. Since the 
accretion disk is unlikely to vary coherently, local fluctuations may be significantly larger than the 
global rms variability. We investigate toy models of quasar accretion disks consisting of a number of 
regions, n, whose temperatures vary independently with an amplitude of gt in dex. Models with large 
fluctuations (gt = 0.35 — 0.50) in 10 2 ~ 3 independently fluctuating zones for every factor of two in ra- 
dius can explain the observed discrepancy between thin accretion disk sizes inferred from microlensing 
events and optical luminosity while matching the observed optical variability. For the same range of 
gt, inhomogeneous disk spectra provide excellent fits to the EST quasar composite without invoking 
global Compton scattering atmospheres to explain the high levels of observed UV emission. Simulated 
microlensing light curves for the Einstein cross from our time-varying toy models are well fit using 
a time-steady power-law temperature disk, and produce magnification light curves that are consis- 
tent with current microlensing observations. Deviations due to the inhomogeneous, time-dependent 
disk structure should occur above the 1% level in the light curves, detectable in future microlensing 
observations with millimag sensitivity. 

Subject headings: accretion, accretion disks — black hole physics — galaxies: active — gravitational 
lensing: micro 



1. INTRODUCTION 

Optical emission from active galactic nuclei (AGN) 
is thought to be due to thermal emission from a 
standard thin accretion disk (|Shakura fe Sunvaevl 119731 : 
iNovikov fe Thornelll973l ). This model is successful in ex- 
plaining many observations of black hole accretion disks, 
such as the peak frequencies of the thermal spectra of 
X-ray binaries and AGN. It is also used to fit contin- 
uum spectra in the X-r ay to measure black hole spin 
(e.g., iShafee et al.l I2006T) . and to infer accretion disk 
size from gr avitational microlen sing; both in absolute 
terms (e.g... IMorgan et all 120101) and as a function of 
wavelength (lAnguita et al.l 120081: IPoindexter et al.ll200l 
lEigenbrod et alj|2008t IMorgan et al.H2010D . Attempts to 
fit continuum AGN spect ra with thin disk models have 
met with mixed success ([Blaes et al.l 120011 : iDavis et"aLl 
12007ft . While multi-temperature blackbody spectra gen- 
erally fit well in the optical, at UV wavelen gths quasar 
spect ra are brighter than predicted (e.g., Bla es et ail 
2001). The X -ray size ha s been shown by m i crolensing 
to be extremely compact (jChartas et al.l l2009: Da i et all 
I2010T ) , disfavoring spectral models invoking global Comp- 
tonizing atmospheres to explain the UV emission. 

Quasar optical variability is also difficult to explain in 
the context of the thin disk model. Long term moni- 
toring of AGN has found almost simultaneous variabil- 
ity ac ross optical wavelengths, with lags of less than 1-2 
davs dCutri et al.lll985t iClayel et alJll991t iKorista et al.l 
fl99l IWanders et al.l Il997t ICollier et"a!Tll998D . Com- 
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paring these lags to the radii dominating the thin disk 
emission at these wavelengths gives a traveling speed of 
0.1c for the variability m echanism (|Krolik et al.l 119911 : 
ICourvoisier fc Cravellll99ll ). This would force the vari- 
ability in AGN to be communicated at the local sound 
or Alfven speed rather than on the infall timescale asso- 
ci ated with disk insta b ilities . 

iWambsganss et ail (Q~990) and iRauch fc~ Blandford 
(1991) first measured the size of an accretion disk from 
microlensing, and found that accretion disks were smaller 
than expected from the thin disk model for the observed 
L„. In recent years, the opposite trend has e merged 
(iPoolev et al.ll2007l: iDai et aLlteoiCt IMorgan et aL|[2010t 
iBlackburne et al.l |2010| )~ Sizes are now robustly found 
to be a factor of ~ 4 larger on average than expected 
from the thin disk model at IR/optical/UV wavelengths. 
The size vs. wavelength relation predicted by the thin 
disk, r cx A 4 / 3 , is within the large range allowe d by mi- 
crolensing o bservations (lEigenbrod et al.ll2008f). At the 
sam e time. iKellv et all (120091 ) . iKo zlowski et all (|2010l ) 
and lMacLeod et aTT i 20101 ) have studied large samples of 
quasar light curves. They find that optical quasar vari- 
ability is well described by a damped random walk which 
returns to a mean value on a typical timescale of 200 days 
with variability amplitudes of ~ 10 — 20%. 

It is unlikely that the observed quasar variability is 
caused by a coherently varying accretion dislQ, but 
rather is probably the added effect of many smaller, inde- 

1 The implications of co herent variations for microlensin g mea- 
surements are discussed by Blackburnc & Kochanck (2010). 
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Figure 1. False color temperature maps of damped random walk models with n = 2200 (left), n = 550 (middle) and n = 140 (right). 
The images are interpolated onto a Cartesian grid for plotting, and the color scale increases logarithmically from blue to red to yellow to 
white with a dynamic range of 100. 




0.01 0.10 1.00 

o Li / <L V > 



Figure 2. Median size increase vs. relative variability for various 
zone sizes for the damped random walk model. Each sequence has 
values of ctt from 0.1 — 0.8. The allowed region from microlensing 
and variability studies is shaded. The lines show analytic fits to 
the variability from the damped random walk model (Eq. [TJ, and 
to the fractional increase in half-light radius from the log-normal 
model (Eq. |4j. For n > 100 and cry ~ 0.4, damped random walk 
disks can be large enough to explain microlensing while matching 
the observed variability. 

pendently varying regions. Many such model s have been 
propo sed to explain quasar variability (e.g. iLvubarskiil 
1997). In this Letter, we demonstrate that for the ob- 
served variability characteristics, such an inhomogeneous 
disk can simultaneously explain multiple discrepancies 
between AGN observations and accretion disk theory. 
Inhomogeneous disks can be large enough to explain the 
microlensing observations, while their temperature fluc- 
tuations on small spatial scales naturally explain the 
observed simultaneous variability across optical wave- 
lengths. Temperatures exceeding the local thin disk 
value lead to broader spectra extending into the UV, con- 
sistent with quasar spectra without invoking a Compton 
scattering medium. 

2. INHOMOGENEOUS ACCRETION DISKS 

We assume that i) the optical/UV emission observed 
from AGN originates in an optically thick accretion disk, 
ii) Variations in the disk occur locally, and are uncor- 
rected on large spatial scales, hi) Fluid in the disk is 



on circular Keplerian orbits. Assumption i) allows us to 
model the disk emission using temperature alone, while 

ii) requires that the disk be inhomogeneous. Assumption 

iii) is only used for computing microlensing light curves, 
where the disk surface brightness map enters into the 
magnification light c urve produced. 

iKellv et all ()2009D found that quasar light curves are 
well described as a CAR(l) process, a random walk that 
tends to return to a mean value on a typical timescale. 
The timescale they found was roughly 200 days, con- 
sistent with the thermal timescale in the predicted op- 
tical emission region of many AGN. Thi s behavior was 
confir med for SDSS Stripe 82 quasars by [MacLeod et al] 
(2010). The amplitude of variations in typical sources 
is 10 — 20%. To produce this variability amplitude with 
multiple, independent regions in the disk requires larger 
local variations (total variance cx TV -1 for TV indepen- 
dent zones). The local accretion disk flux and effective 
temperature will then no longer be a monotonic function 
of radius r, but will fluctuate with azimuth and time. 
This will cause the disk spectrum within an annulus to 
have stronger emission at shorter wavelengths than if the 
same flux were emitted at constant temperature. Con- 
sequently, the outer portions of the disk will contribute 
more flux than for a uniform disk, causing the disk to 
appear larger at a particular wavelength. 

The thin disk modeQ is broadly consistent with many 
observed properties of black hole accretion flows, and 
follows from the conservation of angular momentum and 
energy in the gravitational potential of the black hole. 
Its global time-averaged properties are likely correct. As- 
suming a geometrically thin, optically thick disk the local 
flux can be written F v (r, 4>, t) = ttB v (T), where F v is the 
flux at the frequency v, B v is the Planck function and 
T = (F/db) 1 / 4 , and Ob is the Stefan-Boltzmann constant. 
In the thin disk case, T cx r~ 3 / 4 well outside the inner 
disk edge and is independent of <fi and t everywhere. 

3. DAMPED RANDOM WALK DISK MODEL 

We allow T to be a function of (j> and t as well as r by 
dividing the disk into n evenly spaced zones in log r and 
<f> per octave in radius. Motivated by variability stud- 

2 In this Letter, we us e "thin disk model" to mean that of 
Shakura & Sunvacv (1973). Inhomogeneous disks may also be ge- 
ometrically thin. 
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Figure 3. Top: Sample microlcnsing light curves for all mod- 
els, shifted for legibility. The mean of each curve is ~ 1.3. Best 
fit power-law disk model light curves are shown as dotted lines. 
Bottom: Fractional difference between inhomogeneous and best fit 
power-law disk magnification light curves. Deviations are largest 
at small n and large or. 

ies, we let log 10 T follow an independent damped ran- 
dom walk in each zone with an amplitude <jt that is 
independent of radius. The local temperature returns to 
a mean value on a characteristic timescale of 200 days 
that is also taken to be independent of radius, and the 
mean value is chosen so that (cr&T 4 ),^ = F(r), where 
F(r) is the thin disk flux at the radius r. The results 
in this paper are completely insensitive to the choice of 
characteristic timescale and its radial dependence. The 
total flux is calculated by integrating over the area of 
the disk assuming face-on viewing and ignoring all rela- 
tivistic effects. The accretion flow undergoes Keplerian 
rotation, but the zones remain completely independent 
at all times. Sample temperature maps from this model 
for a few values of n are shown in Figure [1] 

Figure [2] shows the ratio of half-light radii between the 
damped random walk disk and a standard thin disk at 
the same flux as a function of the variability, <jl v - The 
overall variability is fit well by the expression, 



8n y 



(1) 



For a large number of zones, the disk can be signif- 
icantly larger than a thin disk while producing little 
variability. For few (large) zones, the overall variability 
becomes unrealistically large before the accretion disk 
becomes large enough to explain the microlensing mea- 
surements. By design, power spectra and structure func- 
tions from this damped random walk model are consis- 
tent with studies of quasar variability dKellv et al.ll2009t 
iKozlowski etalll20lol: iMacLeod et al.ll2010j ). This model 
is by no means unique. We have also tried other time- 



variable toy models for inhomogeneous disks, which can 
reproduce the observed variability characteristics. All 
models considered give similar size- variability amplitude 
relations. 

In the limit n — > 00, the damped random walk model 
becomes time-independent with a log-normal distribu- 
tion of disk temperatures in each annulus whose variance 
is (In 10 ctt) 2 /2. In this limit, the flux from an annulus 
can be written, 
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where z = hv/kT, T is the thin disk temperature, and 
the local temperature in an infinitesimal piece of the an- 
nulus is wT. The specific luminosity can be computed 
from, 



2-ktF v (t\ ax)dr. 



(3) 



The size increase as a function of ut is the ratio of 
half-light radii from this model assuming T cx r~ 3 ' 4 to 
those computed from L„(0). Sizes are compared at fixed 
L v since the half-light radius, rt, is a function of L v 
(r/j oc ignoring the disk inner edge). The relation, 



rh(cr T )/r h (0) 
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gives an excellent fit to numerical calculations. In models 
with or = 0.35 — 0.5, 5 — 10% of the zones in a given 
annulus produce 50% of its flux in agreement with the 
analytical log-normal (n — >■ 00) model. The fractional 
area contributing half the flux decreases with increasing 
n and ctt- 

3.1. Microlensing Simulations 

Thus far, we have compared ratios of inhomogeneous 
to smooth model thin disk half-light radii at the same 
value of L v . However, microlensing observations do not 
measure half-light radii. Instead, they measure uncor- 
rected variability between different images of strongly 
lensed quasars, caused by the motion of stars in the lens 
galaxy. Sizes can be measured by assuming a power-law 
thin disk with T(r) cx r~ 3 / 4 and fitting magnification 
light curves with surface brightness m aps parameterize d 
by the radius, r s , where hv — kT{r s ) (|Kochanek|[2004). 
For the power-law disk, r s can also be calculated from 
the observed L v in the absence of all lensing effects. The 
discrepancy between the values of r s inferred from mi- 
crolensing versus r s inferred from L v is 0.6 ± 0.3 dex 
([Morgan et al.l20Tlil l. In the power-law disk, = 2.44r s . 
We would like to determine whether microlensing light 
curves from our inhomogeneous toy model can be fit well 
by the power-law disk model, and whether the values of 
r s inferred from microlensing simulations of the inhomo- 
geneous model are consistent with our previous results 
for the half-light radius. 

Using the code described by [Wambsganss (1999|), we 
simulate magnification patterns for QSO 2237+0305, the 
Einstein Cross, u sing lens galaxy parameters found by 
iKochanekl (|2004f) . Choosing random orientations and 
starting positions, many magnification light curves are 
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Figur e 4. Top panel: M odel spectra compared to the composite 
from Zheng ct aQ {1993). The open diamonds show a thin disk 
spectrum. Bottom panel: Model spectra (lines) compared to thin 
disk spectrum (diamonds) for an annulus at r = 20r ms , where 
r ms is the inner disk edge. The black hole mass is 10 8 Mq with 
L/L edd = 0.1. 

produced from the time-varying toy models at multiple 
observed frequencies, and for the power-law disk at many 
values of r s . The microlensing light curves are then fit 
to determine best fit values of r s . Sample light curves 
are shown in the top panel of Figure [3] along with best 
fit power-law disk light curves, while the bottom panel 
shows fractional deviations between the two. 

The time-steady, power-law disk model produces excel- 
lent fits to the microlensing light curves. In the damped 
random walk model, r s from microlensing is consistent 
with the half light radius of the disk at various frequen- 
cies by a fixed fraction that depends on n. For many 
(few) zones, r s underestimates (overestimates) the time- 
averaged disk half-light radius. The largest deviations 
are of order 30% and will not dominate the error in esti- 
mating quasar accretion disk sizes from microlensing. 

3.2. Spectra and Light Curves 

Damped random walk model spec tra are compared to 
the composite quasar spectrum from lZheng et al.l (1997) 
in the top panel of Figure |4j The temperature normal- 
izations have been chosen to provide the best fit. In all 
cases, the toy model spectra have larger F\ at smaller A 
than the thin accretion disk model. The bottom panel 
of Figure S] shows local spectra at r = 20r ms , where r ms 
is the inner disk edge. The local inhomogeneous spectra 
peak at a smaller value of XF\ than the single tempera- 
ture blackbody. To produce the same total flux at that 




Figure 5. Allowed 68% regions of the o~t vs. n parameter space 
from microlensing (blue), variability (red) and spectral (orange) 
observations. Their intersection is shaded green. Contours of 1% 
(solid) and 2% (dotted) rms deviation between model microlensing 
light curves and best fit power-law temperature disks are overplot- 
ted. 

wavelength, the emission must arise from a larger area. 
Thus, the disk appears larger at that wavelength. 

These spectral differences will cause parameter values 
inferred from optical spectra to differ between the two 
models. We find that inferred black hole masses (accre- 
tion rates) are smaller (larger) when fitting inhomoge- 
neous spectra with thin disk models. The discrepancy is 
typically less than a factor of two, and is unlikely to be 
a significant source of error in estimating accretion disk 
parameters from spectral fitting. 

The damped random walk model produces simultane- 
ous optical variability, in that the cross-correlation be- 
tween all wavebands peaks at a lag of zero timesteps. 
The lag is zero regardless of the chosen timestep, because 
neighboring zones are completely independent. Variabil- 
ity is stronger at higher frequencies, in qualitative agree- 
ment with observations. Quantitative predictions for this 
relation, for cross-correlations between light curves at dif- 
ferent observed wavelengths, or for trends with luminos- 
ity or black hole mass will require a physical model for the 
variability mechanism. These quantities are highly sen- 
sitive to the particular inhomogeneous disk prescription 
chosen, many of which satisfy the existing observational 
constraints. 

3.3. Combined Observational Constraints 

Figure [5] summarizes our main results. The 68% al- 
lowed regions of the o~t vs. n parameter space from 
microlensing, spectral and variability measurements are 
shown for the damped random walk model, along with 
contours of rms deviations between model and best fit 
power-law disk microlensing light curves at 2500A. The 
disk size and variability constraints are the same as 
in Figure [2] The spectral constraint is found from 
fitting mode l spec tra to the composite spectrum from 
iZheng et al.l (fl9 97i . varying the temperature normaliza- 
tion in each case to obtain the best fit. 

The independent spectral and microlensing size con- 
straints require similar amplitudes of temperature inho- 
mogeneity. Producing the observed level of variability 
requires larger <jt at larger n. All constraints intersect 
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in the green region, giving best fit parameter values of 
0.35 < a T < 0.5 and 100 < n < 1500. The models pre- 
dict short timescale deviations from smooth power-law 
disk models in microlensing light curves at or above the 
1% level on average, with maximum deviations of 10% 
(see Figure [3]). These are likely to be lower limits. De- 
viations from other toy models tried can be a factor of a 
few larger. 

4. PHYSICAL MECHANISMS 

Quasar accretion disks are expected to be inhomoge- 
neous to some degree. However, to produce the ob- 
served sizes and broad spectra observed in AGN, the 
inhomogeneities must be large - factor of 3 — 5 vari- 
ations in temperature at a given radius. Disk insta- 
bilities are the most promising mechanism for produc- 
ing inhomoge neous disks. The magnet orotational insta- 
bility (MRI, IBalbus fc Hawlevl 119911) is now believed 
to operate in a wide range of accretion flows, caus- 
ing the disk to become fully turbulent and providing 
a means for outward a ngular momentum tr ansport and 
accretion. Local (e.g., IMiller fc' Stone 20Q3) and global 
simulations of thin MRI disks have been performed 
both in pseudo- Newtonian poten tials (jArmitage et al.l 
120011: I Armitage fc Reynolds! 120031) and recently i n full 
gener al relativity (iShafee et al.l l2008t iNoble et al.l 120091 
l2010t iPenna et al.l 120101 ) . None of these simulations in- 
cluded radiation forces, which likely dominate the dy- 
nami cs of AGN accretion disks ov er a large ran g e in ra- 
dius (jShakura fc Sunvaevl Il97l . IHirose et all {2009a) 
studied a set of local simulations in radiative MHD, 
many of which were radiation-dominated. Their sim- 
ulations show temperature variations with factors of 
~ 2 in a local patch of disk . Global general rela- 
tivistic MHD simulation s from iFragile et "aLl (|2007l ) and 
iMcKinnev fc Blandfordl (pool ) have gas temperatures 
that can vary by factors of 2 — 3 at a single radius in the 
midplane, but are considerably thicker than AGN disks 
and neglect radiation altogether. These distributions 
correspond to <tt — 0.1 — 0.2. It is unclear whether the 
MRI-driven fluctuations in radiation-dominated global 
disks will be larger. 

Thin, radiation-dominated disks have long been be- 
lieved to be subject to a strong thermal instability, 
since the heating rate depends on a hig her power of 
the tem perature than the cooling rate. IHirose et al.l 
(2009b) find that the dynamics of shearing box simu- 
lations of radiation-dominated disks are well described 
by the alpha model on long timescales. The flow 
is highly variable, but thermally stable. Radiation- 
domin ated disks are also subjec t to an inflow insta- 
bility (jLightman fc Eardlevl 11974 ). Radially extended, 
radiation-dominated MHD simulations are needed to de- 
termine whether this instability operates in AGN. If 
so, it would likely create large local temperature gra- 
dients in the accretion flow between hot, tenuous gas 
and cooler optically thick regions. Finally, it is possi- 
ble that a photon bubble instability could be present i n 
radiation-dominated accretion disks ([Turner et al.l l2005). 
which can cause the flux escaping the disk to rise by 
factors of several locally. This may prevent inhomoge- 
neous disks from becoming dynamically unstable, despit e 
exceeding the Eddington limit locally (Begelman 2006). 
The disk would also remain stable if magnetic pressure 



prov ides the vertical support i n quasar accretion disks 
(e.g. [Begelman & Pringlc 2 0~07l) . Alternatively, the disk 
may drive a wind when the Eddington limit is exceeded. 

5. CONCLUSIONS 

Standard thin disk accretion has formed the basis for 
understanding X-ray binaries and AGN for nearly 40 
years. However, it has always been difficult to extend 
the model to explain optical quasar variability. In re- 
cent years, microlensing observations of multiply im- 
aged quasars have provided a probe of the disk struc- 
ture, finding that quasar microlensing sizes are robustly 
larger than the flux sizes predicted from thin disk the- 
ory. If the average temperature structure remains iden- 
tical to that in thin disk theory but is highly inhomoge- 
neous, accretion disks can be large enough to explain the 
sizes found by microlensing while matching the observed 
level of optical variability. The level of inhomogeneity 
(cr T = 0.35 — 0.50 with n — 10 2 ~ 3 ) required to explain 
the discrepancy in microlensing sizes is in excellent agree- 
ment with that necessary to fit observed quasar spectra. 
Such inhomogeneous structure produces short timescale 
variations in microlensing light curves that should be 
larger than ~ 1%. The amplitude of the temperature 
fluctuations can be further constrained from measuring 
the size of these deviations. The range of temperatures 
in small regions of the inhomogeneous disk explains the 
simultaneous variability observed across optical wave- 
lengths. 

We have demonstrated this idea with an unphysical toy 
model. Proper modeling of an inhomogeneous disk will 
require global MHD simulations of radiation-dominated 
accretion disks. It is unclear whether the MRI alone is 
sufficient to produce the required temperature fluctua- 
tions, or whether additional disk instabilities or other 
variability mechanisms are also important. 
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